Discovery of new quasi-periodic oscillations in the X-ray transient 

source V 0332+53 
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ABSTRACT 

We report the discovery of a new quasi-period oscillation (QPO) at 0.22 Hz, centered on the 
source spin frequency of the high mass X-ray binary system V 0332+53 when the source was 
observed during its November 2004/March 2005 outburst by RXTE. Besides this new QPO, we 
also detected the known 0.05 Hz QPO. Both the 0.22 and 0.05 Hz QPOs stand out clearly at a 
mid-flux level of the outburst within January 15-19 2005, and later at an even lower flux level as 
the width of 0.22 Hz QPO drops. No evolution of the centroid frequency with the flux is seen in 
either QPO. The rms value below 10 keV is around 4-6% for both QPOs and decreases at higher 
energies. We discuss our results in the context of current QPO models. 

Subject headings: binaries: general — stars: individual (V 0332+53) — X-rays: binaries — X-rays: 
individual (V 0332+53) — X-rays: stars 



1. INTRODUCTION 

The high mass X-ray binary (HMXB) tran- 
sient V 0332+53 (X 0331+53) was discovered dur- 
ing a bright outburst phase in July 1973 with 
the Vela 5B observatory (Terrell & Priedhorsky 
1984). After the discovery, the source was ex- 
tensively studied in three new outbursts: in 1983 
by Tanaka et al. (1983) with EXOSAT, in 1989 
by Makishima et al. (1990a) with Ginga and the 
last one in 2004/2005 by Zhang et al. (2005) and 
Kreykenbohm et al. (2005) with RXTE and IN- 
TEGRAL. No outburst was detected by BATSE 
on board CGRO during its era between 1991 and 
2000 nor by RXTE until November 2004 when the 
source outburst was caught by the all sky monitor 
(ASM). 

The recurrent HMXB transient system consists 
of an accreting neutron star (NS) and an early- 
type 08-9 Ve optical counterpart BQ Cam lo- 
cated at a distance of ~ 7 kpc (Bernacca et al. 
1984; Honeycutt & Schlegel 1985; Negueruela et 



al. 1999). The pulse period was flrst measured at 
4.375 s with a 34 day orbit period, and a moder- 
ately eccentric orbit e 0.3 (Stella et al. 1985). 
The morphology of the pulse profile changes as 
a function of energy. It shows a prominent and 
broad single-peak at lower-energy and a secondary 
peak which becomes more evident above 6 keV 
(Unger et al. 1992; Kreykenbohm et al. 2005). 
However, Stella et al. (1985) and Zhang et al. 
(2005) found that the pulse profile is also depen- 
dent on the source luminosity - evolving from a 
double pulse when the source is bright to a single 
pulse when the source fiux is lower. The period 
was found to decrease during outburst, with a re- 
cently measured value of —P ^ 9 x 10~^^ s s~^ 
(Zhang et al. 2005). Possible evidence of an accre- 
tion disk was found with the detection of a 0.05 
Hz QPO (Takeshima et al. 1994). 

Observations (Makishima et al. 1990b; Coburn 
et al. 2005; Soldi et al. 2005; Kreykenbohm et al. 
2005) show that the spectrum is best fitted by an 
absorbed power-law with an exponential cut-off 
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and the additional components of cyclotron ab- 
sorption lines, located at 25 keV, 50 keV, and 70 
keV. The magnetic field at the NS surface is esti- 
mated to be around 2.7x10^^ Gauss. 

In this letter we report the discovery of a 
new QPO based on the RXTE observations of 
V 0332+53 during the 2004/2005 outburst. 

2. OBSERVATIONS and DATA 

After the initial follow-up observation by RXTE 
starting November 27, 2004, Target of Opportu- 
nity (ToO) observations were performed between 
December 28, 2004 and February 15, 2005, the 
data from which were made public. The RXTE 
observations caught the source V 0332+53 at the 
outburst maximum on December 28-30, 2004 and 
during the outburst decay phase on January 4- 
6, 15-19, and February 12 15, 2005. In Fig.l we 
show the four time intervals in which the RXTE 
observation was performed. Hereafter we call the 
four datasets /, //, /// and IV, respectively. For 
our timing analysis, we used these datasets from 
the Proportional Counter Array (PGA; 2-60 keV) 
(Jahoda et al. 1996). The PGA data were col- 
lected in the binned and event mode, where the 
default selection criteria were applied. We used 
only the PGU2 which was on during the whole 
observation; the net PGU2 exposure time was 155 
ks. The data reduction was performed using the 
standard software package FTOOLS version 5.3.1, 
and the background contribution was subtracted 
using the pcabackest version 4.0. 

3. RESULTS 

3.1. ASM AND PCA LIGHT CURVE 

In Fig.l we show the one day averaged ASM 
light curve of the entire outburst in the 1.5-12 
keV energy band. We converted the ASM count 
rate to flux using 1 Grab « 75.5 cts/s (Remillard 
2004). The true flux is certain to be substantially 
higher since the source is extremely hard, even for 
X-ray pulsars. The outburst started on Novem- 
ber 23, 2004 (MJD 53332) and reached the maxi- 
mum flux level of 1.2 Grab on December 29, 2004 
(MJD 53368.6+0.4). The flux then decayed with 
an exponential time-scale, F oc e~*/^^-^ , to quies- 
cence. The first evidence of a new X-ray outburst 
of this source could be traced back to January 2004 
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Fig. 1. — Top panel: one day averaged 
RXTE/ ASM light curve in the 1.5-12 keV en- 
ergy band of V 0332+53 during the 2004/2005 
outburst. The vertical arrows indicate the times 
of the RXTE/PCA observations. Middle panel: 
hard color (12.0-5.0) keV to (3.0-5.0) keV. Bot- 
tom panel: soft color (3.0-5.0) keV to (1.5-3.0) 
keV. 

through observed optical brightening (Goranskij & 

Barsukova 2005). To analyze the V 0332+53 spec- 
tral variability as function of the outburst time, we 
defined two colors. The soft color is the ratio of 
the 3.0-5.0 keV band to the 1.5-3.0 keV band and 
hard color is the ratio of the 5.0-12.0 keV band 
to the 3.0-5.0 keV band. Golors are built from 
the one day averaged ASM light curve. From the 
color ratios it is evident that the source under- 
went a spectral change during the outburst. We 
find that the hardness in the soft color increases 
to a maximum at the outburst peak, and that the 
transition in hardness occurs at an earlier time for 
the hard color. 

In Fig. 2 we show the 0.5 day averaged PGA in- 
tensity evolution of the source in the 2.0-16.5 keV 
energy band. Using the 16 s time-resolution Stan- 
dard 2 mode PGA data we calculated the color 
ratios. The PGA soft color (3.2-5.7 keV/0.7-3.2 
keV) shows a linear trend consistent with the ASM 
one in the same time interval. The PGA hard color 
(9.4-16.5 keV/5. 7-9.4 keV) shows a spectral hard- 
ening during the outburst flux decay, compared to 
the little change of the spectral shape in the soft 
color. 
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Fig. 2.— The PCA light curve (2.0-16.5 keV) (top 
panel) and the time evolutions of soft color (mid- 
dle panel) and hard color (bottom panel). Here 
MJD 53367 corresponds to December 28, 2004. 

3.2. TIMING ANALYSIS 

We searched for QPOs in the 2-60 keV energy 
band, using binned mode and event mode PCA 
data. We computed separately for each of the 
four datasets a Power Density Spectrum (PDS) 
in the frequency range between 1/256 and 8 Hz 
from Fast Fourier Transforms. The PDS is nor- 
malized using the method described by Miyamoto 
et al. (1991), and the Poisson white noise is sub- 
tracted. In the resulting PDS two evident QPO 
signals are present at ^0.05 and ^^0.22 Hz. The 
QPO feature at 0.22 Hz is newly detected; the 
low frequency QPO at 0.05 Hz was previously ob- 
served with Ginga (Takcshima ct al. 1994). Both 
QPOs are detected in the dataset ///at very high 
significance and in the dataset IV at lower statis- 
tical significance. Outside these two time periods 
the QPOs are not obvious. The PDS spectrum 
also shows high significant peaks at the spin fre- 
quency, corresponding to a fundamental frequency 
of 0.228 Hz, and its higher harmonics. The new 
QPO is very close to the spin frequency of the NS. 
This provides the first detection of a clear QPO 
riding on the spin frequency of a neutron star for 
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Fig. 3. — The power density spectra of V 0332-1-53 
obtained during the dataset 777 and IV. The spikes 
were removed when the PDS was fitted (see text). 

The accurate centroid frequency and the 
quality factor Q = uo/dv of the QPOs were found 
by fitting the PDS with two Lorcntzian functions. 
The broad band-limited noise component was fit- 
ted using a power-law and a Lorentzian func- 
tion. The best fit was also found by removing 
the spikes of the spin frequency and its harmon- 
ics (3-5 frequency bins for each spike) from the 
PDS. For dataset /// (January 15-16, 2005), we 
found the centroid frequency for the new QPO 
at uo = 0.223 ± 0.002 Hz and the quality fac- 
tor Q w 2.1. The lower frequency QPO was 
at uo = 0.049 ± 0.007 Hz with a quality factor 
of Q » 1.3. One month later, using the whole 
dataset IV we found for the higher frequency QPO 
I/O = 0.227 ± 0.002 Hz and Q ^ 12. 

The Fourier transform of the time series re- 
sulted in the PDS with best fits shown in Fig. 3. 
In Fig. 4 we show the evolution of the QPO cen- 
troid frequencies and widths versus the X-ray flux 
for the datasets /// and IV. We found that the 
centroid frequencies of both QPOs do not change 
significantly with the flux, while the width of the 
0.22 Hz QPO drops with lower flux level. 
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Fig. 4.— Flux evolutions of the centroid frequen- 
cies and widths for the 0.22 Hz QPO (first two 
panels) and for the 0.05 Hz QPO (lower two pan- 
els). 

To investigate the QPO energy dependence, we 
estimated the fractional root-mean-square (rms) 
from the PDS. We consider only the dataset /// 
for this analysis, since the two QPOs stand out 
significantly in this outburst phase. We find for 
both QPOs that the rms is roughly constant below 
10 keV and decreases by a factor of ~2 at higher 
energies (Fig. 5). The fact that this trend was not 
found in the 1989 outburst using Ginga data may 
be due to the relatively small effective area of the 
detector (Takeshima et al. 1994). 

4. DISCUSSION 

We have discovered a new QPO in V 0332+53 
at vq = 0.223 ± 0.002 Hz, close to the spin fre- 
quency of the neutron star. This new QPO was 
visible in the PDS at the same time as the well 
known i/q = 0.049 ± 0.007 Hz QPO. The fact that 
the QPO is found near the spin frequency sug- 
gests that it is produced on or near the surface 
of the NS. QPOs found in HMXB pulsars have 
a frequency range of about 10 mHz to 400 mHz 
(Finger 1998; Ghosh 1998, for review ) , though 
none have so far been detected at the NS spin fre- 
quency. Broadening features in the PDS at the 
spin frequency were detected as well in some other 
HMXB such as Vela X-1, although they were not 
classified as QPOs (Lazzati & Stella 1997). The 
QPOs found in HMXBs are considerably lower 
compared with those found in low mass X-ray bi- 
nary (LMXB) systems which lie in a frequency 
range from ~ 6 to 50 Hz for the low frequency 
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Fig. 5. — The energy evolutions of the rms for the 
0.22 Hz QPO (upper panel) and for the 0.05 Hz 
QPO (lower panel). 

QPO and - 300 - 1200 Hz for the higher fre- 
quency QPO (e.g., van der Klis 2004, for review). 
In LMXB, one of the QPOs is interpreted as the 
NS spin frequency; this may imply a similar ori- 
gin for the newly observed QPO in the HMXB 
V 0332+53. 

The precise nature of the QPO is still unknown. 
However, a series of models have been proposed 
to explain the QPO features observed in compact 
X-ray binary systems (e.g., van der Klis 2004; 
Finger 1998; Ghosh 1998, for review). Among the 
most frequently used are the magnetospheric beat 
frequency model (MBFM) and the Keplerian fre- 
quency model (KFM). In the MBFM, the QPO 
frequency is interpreted as the beat frequency be- 
tween the stellar rotation, Vg, and the Keplerian 
rotation, vk, at the inner edge of the accretion 
disk, i.e., vqpo — i^K — t^s (Alpar & Shaham 1985; 
Shibazaki & Lamb 1987). The X-ray flux is pro- 
portional to the mass accretion rate, i.e., F oc M. 
In the KFM (van der Klis ct al. 1987), the QPOs 
arise from the modulation of the X-rays by in- 
homogeneities in the inner disk at the Keplerian 
frequency, i.e., vqpo = vk, and the X-ray flux is 
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Both models predict a systematic change of the 
QPO frequency with the mass-accretion rate or 
with the flux. However, as shown in Fig. 4, such 
a feature is not observed in V 0332+53. From 
dataset III to IV the X-ray flux changed by a fac- 
tor of ^ 4.5, but the centroid frequencies of both 
QPOs remain constant within the error bars. This 
suggests that neither the MBFM nor the KFM 
can explain the QPO phenomena as observed in 
V 0332+53. Similar conclusions have been drawn 
by different authors. For example, the expected 
luminosity related to the QPO frequency does not 
always agree with the observed luminosity (Fin- 
ger 1998). Other observations have not been sup- 
ported by the MBFM model, e.g., as found by 
Takeshima et al. (1994) through the comparison 
of the estimated Alfven radius from accretion disk 
theory (Ghosh & Lamb 1979a,b). 

One alternative model to explaining such cor- 
relations was given by Lazzati & Stella (1997) and 
Burderi et al. (1997). They argued that the ac- 
cretion flow is inhomogeneous near the surface of 
the neutron star, and random shots are character- 
ized by an arbitrary degree of modulation. The 
combination of rotation and radiative transfer ef- 
fects should therefore produce a periodic modula- 
tion of the shots similar to that of any continuum 
X-ray emission from the polar caps. They pro- 
posed that a coupling between the periodic and 
red-noise variability should be frequently present 
in X-ray pulsars. The discovery of a new QPO 
with the frequency of the spinning neutron star 
provides an additional example which could be un- 
derstood within the scenario of this model. 
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ence Archive Research Center Online Service, pro- 
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